Abstract. We investigate the effect of pair creation on a 
shock structure. Particles, accelerated in the shock via the 
first order Fermi process, are supposed to cool by inverse 
Compton process on external soft photons, resulting in a 
cut-off power law shape of the particle distribution func- 
tion. The high energy photons produced are thus able to 
create pairs, through photon-photon annihilation. The in- 
crease of the pair pressure may then be sufficient to modify 
the shock profile. We show that there is even a limit of the 
t-H ' pair pressure (of the order of 20% of the ram pressure of 

. the upstream flow) above which the shock cannot exist 

| any longer. 

. Conversely, significant changes of the flow velocity pro- 
file will also modify the spectral index and the high en- 
ergy cut-off of the particle distribution function. Hence 



the number of particles able to trigger the pair creation 
, process will change, modifying the pair creation rate ac- 
t-H ' cordingly. Taking into account these different processes, 
we self-consistently derive the flow velocity profile and 
' the particle distribution function. We show that, in some 
, region of the parameter space, the system can converge 
' towards stationary states where pair creation and hydro- 
y—i , dynamical effects balance. 

lO ' We discuss the application of this model to explain the 
, high energy emission observed in compact objects. We 
! show that hard X-ray spectra (ax < 1.) are only obtained 
--.^ for small pair pressure and we don't expect any strong 
i-£h annihilation line in this case. We suggest also a possible 
variability mechanism if the soft photon compactness de- 



Or 
i 



O ' pends itself on the pair density of the hot plasma, such as 



in 
c3 



expected in reillumination models. 
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1. Introduction 

It is well-known that particles can be accelerated up to 
very high energies by crossing a magnetized shock front 
many times and thus experiencing the so-called first 
order Fermi process (Axford et al., 1976; Krimsky 1977; 
Bell p.978} Blandford & Ostriker (1978Q. The simple linear 



theory may be however insufficient to describe a realistic 
situation because the accelerated particles may induce 
a strong non-linear back reaction on the shock itself. 
For example they may produce magnetohydrodynamic 
(Jokipii 1976f Lacombe 1977) or electrostatic (Nishikawa 



et al. 1994) instabilities, and also strongly alter the 



hydrodynamics of the shock through their own pressure 



(Drury et al. 1982; Schneider & Kirk 1987) 



In the case of astrophysical objects like the Active 
Galactic Nuclei (herafter AGNs), where shocks are 
particularly attractive processes to produce in situ high 
energy particles, the presence of important radiation 
fields requires also to take into account particles-photons 
interactions in the shock region. A large number of works 
have already studied particles acceleration at shocks 



including radiative cooling processes (Webb et al. 1984 
Schlickeiser 11984; Bregman 1985; Heavens & Meisen- 



heimer 1987; Biermann & Strittmater 1987; Protheroe 



& Stanev |1999| ; Drury et al. |F399|) or the feedback of 
the accelerated particle pressures on the shock structure 
(Blandford |1980| ; Drury & Volk |1981| ; Heavens |1983| ). 



The aim of this paper is to investigate also the effect 
of pair creation on the shock structure. Indeed particles 
accelerated by the shock can be sufficiently energetic to 
boost, via Inverse Compton (hereafter IC) process for 
example, surrounding external soft photons above the 
rest mass electron energy and thus allowing the trigger of 
the pair creation process. The increase of the associated 
pair pressure is thus able to disrupt the plasma flow and 
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eventually, for too high a pressure (due to catastrophic 
run-away pair production for example) to smooth it 
completely. Reversely, significant changes of the flow 
velocity profile may modify the distribution function of 
the accelerated particles in such a way that the number 
of particles enable to trigger the pair creation process 
can change, modifying consequently the pair creation rate. 

This non-linearity of the problem is in addition to 
the intrinsic non linearity of the pair creation process, 
which may be more or less important depending on the 
compactness of the soft and hard photons field in the 



source (Svensson 1987). The pair creation process is also 



naturally non-local in the sense that pairs are produced 
by photons which are themselves emitted by particles 
located in an other part of the plasma. The complete 
3D-treatment of this problem is out of the scope of this 
paper. We propose however to use simplifying assump- 
tions in ID-geometry, to study the global behavior of a 
shock structure in an environment dominated by pairs. 

This paper is divided as follows. We first present in 
section 2 the basic hypotheses we use to treat (relatively 
easily) the problem without loss of generality. Some im- 
portant definitions and descriptions of the geometry of 
our model are detailed in section 3. We then present the 
main equations describing the behavior of a shock in a 
pair environment in section 4. In section 5, we first dis- 
cuss the case without hydrodynamical feedbacks on the 
pair process, deducing a limit of the pair creation rate 
above which the shock disappears. When feedbacks are 
taken into account, we study in section 6 the existence of 
stationary states where pair creation and hydrodynamical 
effects balance. We discuss the results in section 7 before 
concluding. 

Throughout the paper upstream and downstream quanti- 
ties are denoted by subscripts 1 and 2 respectively. Par- 
ticles and photons energies are expressed (unless pointed 
out) in unit of m e c 2 , speeds in unit of u%, the upstream 
flow velocity, lengths in unit of the upstream diffusion 
length Di/ui, where D\ is the upstream diffusion coef- 
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ficient, and times in unit of D\/u\. Variables in reduced 
units are overbraced by a ~. 



2. Basic hypotheses 

2.1. Simple geometry 

We assume a ID geometry (the x axis being align with 
the flow velocity cf. Fig. |l|) meaning that the different pa- 
rameters characterizing the flow are homogeneous in each 
section perpendicular to the x axis. However, we will use 
the transverse choc radius R to evaluate the size of the 
pair creation region i? 77 in section 3.2.2. The 1-D assump- 
tion is certainly well justified in the central part of the 
flow where boundary effects are negligible. We also make 
the assumption of a parallel shock, i.e. the magnetic field 
lines are supposed to be parallel to the normal of the shock 
front. We however assume the magnetic field to be slightly 
perturbed near the shock and we admit these perturba- 
tions to be dominated by Alfven waves. It is thus possible 
for particles to be scatter ed by these wav es thr ough pitch 
angle scattering (Jokipii 1976 ; Lacombe 1977 ) and thus 
going back and forth across the shock. 

2.2. Constant coefficient diffusion 

Rigorously the upstream and downstream spatial diffu- 
sion coefficient, D\^ xx and D2, XX , depend slightly on 7. For 
small Larmor radii of the particles compared to the large 
scale of the turbulence (w hich is well-verified for leptons) 
one can show (Casse et al. 2001 ) that Di jXX oc 7 2_/3 where 
(3 — 5/3 for Kolmogorov spectrum. The dependence on 7 
being weak, we then consider Di yXX to be independent of 
the Lorentz factor 7 of the particles. We also assume Di iXX 
to be independent of x. Since we work in ID-geometry, 
the upstream and downstream coefficient diffusion will be 
simply noted D\ and D2 respectively in the following. 



2.3. Plasma dominated by the relativistic pressure 

We suppose the plasma pressure P to t to be dominated by 
the pressure P re i of the relativistic leptons population. 
This seems relatively reasonable since the leptons accel- 
eration time scale is expected to be relatively small (and 
largely smaller than the proton acceleration time scale, 
Henri et al. 1999). P re i becomes thus rapidly larger than 



the thermal pressure. In this case, the flow velocity pro- 
file make already a smooth transition between the up and 
downstream region but acceleration still occurs (Drury et 
al. 119821; Schneider & Kirk 119871). Besides, with this as- 



sumption, the specific heat ratio of the accelerated parti- 
cle plasma is equal to 4/3. Consequently, in the case of a 
strong shock, the compression ratio reaches the value of 7. 



2.4- Assuming a large Alfvenic Mach number 

Although the magnetic field is necessary to produce the 
stochastic acceleration, we will assume that the Alfven 
velocity Va — B/^/p (where B is the magnetic field am- 
plitude and p the plasma density) is much smaller than 
the upstream flow velocity u\ i.e. the Alfvenic Mach num- 
ber Ma = ui/Va is very large. This leads to important 
simplifications in the treatment of the acceleration process 



(Vainio & S chlike iser |l998j , |1999| Schlickeiser et al. |1993 
Henri et al. |1999[): 



— the second order Fermi process is negligible compared 
to the first order one since the ratio of the acceleration 
time scale is approximately ti st /t2nd — Ma~ 2 

— The scattering center compression ratio (which con- 
trols the particle distribution spectral index) is close 
to the gaz compression ratio r = U1/U2 

Consequently, in the following, we will neglect the second 
order Fermi process in our treatment and we will only 
use the gaz compression ratio for the computation of the 
particle distribution function (cf. section 3.3). 



2.5. Compton cooling in Thomson regime 

We assume the particles cooling to be dominated by the 
Compton process onto the soft photons, neglecting the 
Synchrotron process. Again it is justified in the case of 
low magnetic energy density. For simplicity, the external 
soft photons field is supposed to be monoenergetic and we 
note e s the corresponding soft photon energy. This is not 
a bad approximation if the soft photons field is emitted by 
a thermal plasma like that probably associated with the 
Blue Bump emission seen in most AGNs. Furthermore we 
will restrict our study to the Thomson regime i.e. e s "f c <1 
where 7 C is the high energy cut -off L orentz factor of the 
particle distribution (cf. section 3.3.3 ). 



3. Some important parameters and definitions 

3.1. The acceleration process 
3.1.1. The acceleration threshold 

Only particles having a Larmor radius r^ comparable to 
a wavelength of the Alfven spectrum will undergo scat- 
terings (Jokipii 1976 ; Lacombe 1977 ) and thus will go 
back and forth many times across the shock. However, in 
e~ —p + plasma as those we deal with, the non-relativistic 
protons limit the Alfven waves spectrum to wavelength 
greater than 2ttVa/lo cp where u cp is the cyclotron pulsa- 
tion of the protons in a magnetic field B. Thus a particle 
of momentum p will be scattered by Alfven waves if it 
verifies: 

p m p V A . . 

rL = W\-TW (1) 

that is p > mpVA- This puts a severe threshold for res- 
onant interactions, especially for the leptons that must 
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Fig. 1. Schematic view of a shock. The shock discontinuity is represented by the vertical bold line. We have also 
indicated the different parameters defining the acceleration and pair creation region (cf. section 3). The acceleration 
region is the region in between the 2 vertical dot-dashed line. The pair creation one is filled in grey. Particles are 
represented by straight arrows and photons by warped ones. Scales are not respected. 



be already very energetic to participate to the diffusion 
processes. Indeed, from Eq. (0) we can deduce the lower 
Lorentz factor for a relativistic lepton to be accelerated in 
a shock: 



^fscat 



m p Va 



(2) 



which could be easily of the order of a few. Several 
processes are known to accelerate electrons up to such 
an energy as the development of a parallel electric field 
component in magnetic reconnections, or waves (magne- 
tosonic or whistler which have much lower particle energy 
threshold for resonant interaction, Ragot & Schlickeiser 
199S). We will assume that such process accelerates the 
particles up to a characteristic Lorentz factor 7 m i„ for 
which the pre-acceleration and cooling times are equal. 
We further assume that "/ m in is above the scattering 



threshold 7 sca t- 

It is worth noting that the high energy power law spec- 
tra observed in compact objects generaly extend down to 
the kcV range. For non-thermal models, it implies that the 
injection threshold is low (7 sca t <10), which is consistent 
with a small Alfvenic velocity 



3.1.2. The acceleration region 

Particles interacting with the shock are those located 
within about one diffusion length of the shock, where the 
upstream and downstream diffusion length are defined by 
(as previously said, i=l, 2 for up and downstream respec- 
tively) : 
A 

Li = — (3) 
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For a velocity profile distort by the presence of pairs, this 
definition should not be rigorously used. However, we ex- 
pect the change of the velocity in the vicinity of the shock, 
due to the pair pressure, to be relatively small, of the or- 
der of the unity (the compression ratio r — U1/U2 can only 
vary between 1 and 7). The region explored by particles 
will thus have a size not very different from that given by 
Eq. (|J) We will call in the following "acceleration" region 
the physical space defined by — L\ < x < +L 2 (cf. Fig. |l|). 

3.1.3. The acceleration time scale 



produced will be able to give a pair of electron-positron if 
at least : 

/ 3 \ 1/2 

e>l i.e. 7> I — J =7^ (5) 

A particle needs thus to have a Lorentz factor higher than 
jth to generate a pair electron/positron. To fix an order of 
magnitude, for a UV bump peaking near 10 eV (a common 
value in AGNs, Walter et al. |1994[ ), Eq. (§) gives ~/ t h - 
300. We will suppose that 7 m i„ < jth so that particles 
need to be accelerated in the shock to be able to trigger 
the pair creation process. 



At each back and forth across the shock front, a relativistic 
particle gain an energy such that: 
dp 4 ui — «2 
p 3 v cos 9\ 

where v is the particle velocity and Q\ is the angle of the 
magnetic field line with respect with the shock normal. 
In our case (a parallel shock) we assume 0\ = 0. 

On the other hand, the residence time t r of this particle 
in the vicinity of the shock is determined by the diffusion 
coefficient: 

If the upstream and downstream diffusion coefficients 
were equal, the downstream contribution would be 
dominant. However, as shown by Vainio & Schlickeiser 
( |l998|) , there is generally an increase of the turbulence 
level in the downstream flow so that D\ is expected to be 
significantly larger than D 2 . We refrain to use too much 
involved calculations and we assume that t r is simply of 

the order of the first term, i.e. t r ~ — 

u{ 

At each crossing, the particle has a probability i] ~ 

— - to escape (Bell, 1978) and for a relativistic particle the 

v 

crossing frequency is roughly v c ~ — . The acceleration 

r\t r 

rate is therefore (with a factor of the order of the unity) : 
(Ap) r - 1 
pAt 3t r 

We can thus deduce the corresponding acceleration time 
scale: 

„£>i 1 Di 



t, 



u\ r — 1 



2 L acc 



(4) 



3.2. The pair creation process 
3.2.1. The pair creation threshold 

A soft photon of energy e s (in reduced units) scattered 
by a lepton of Lorentz factor 7 will be boosted, via IC 

4 2 

process, to an energy e ~ -7 e s . The high energy photon 



3.2.2. The pair creation region 

If the shock region is compact enough and if part of 
its radiation extends beyond 511 keV, the pair creation, 
through photon-photon annihilation, will become impor- 
tant. The pair creation optical depth r 77 depends on pho- 
ton energy e. Photons with e > 2 produced in Compton 
scatterings will produce pairs mostly with a Lorentz fac- 
tor 7 ~ e/2 by colliding with photons with energies of the 
order of 2/e. Consequently, the optical depth for e >2 can 
be written as follows (Zdziarski & Lightman 1985 ): 



3 
R_ 

T, 



(G) 
(7) 



is the typical attenuation length of a 



where n p h(e) is the density of photon with energy e 
per unit volume and dimensionless energy e, and L e — 

photon with an energy e. Since we expect the photon den- 
sity to decrease with energy, L e is a decreasing function 
of e and thus L t < L D where we define: 
L = L t= \ 

On the other hand, the photon density will suffer from 
geometrical dilution after some distance comparable to R, 
the transverse size of the shock (cf. Fig. |l|). Hence, we 
will define the pair creation region as the region included 
between — i? 77 and +i? 77 , where i? 77 is defined by: 

R ^ = irrt (8) 

This definition has the advantage to take into account 
both photon-photon depletion and geometrical dilution ef- 
fects. In the following, we will suppose, for sake of simplic- 
ity, that the pair creation rate is constant in this region. 

3.3. The flow profile and particle distribution 

3.3.1. The "effective" compression ratio 

The definition of the compression ratio like the ratio be- 
tween the far upstream and downstream flow velocity is 
rather unsatisfactory in our case since it will not give a 
real estimate of the velocity change experienced by a rel- 
ativistic particle in the vicinity of the shock, where the 
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particle distribution is principally built. We will thus de- 
fine an "effective" compression ratio, noted simply r, as 
the ratio between the upstream velocity in — lOLi and the 
downstream velocity in +IOL2. We have checked that in 
the case of a strong shock without pairs, this definition still 
gives a compression ratio very near the expected value of 
7. 

3.3.2. The spectral index 

It can be shown that the solution of the evolution equa- 
tion of the particle distribution function including the pair 
creation process (cf. Eq. Jl3|)), still has a energy power law 
dependence, as it is effectively the case without pairs (cf. 
Appendix [A|) . The spectral index s [] keeps also the same 
expression in function of the compression ratio i.e.: 
«(r) = g-2 = (r + 2)/(r-l). (9) 
For plasma dominated by relativistic pressure, the com- 
pression ratio is necessarily smaller than 7 so that s(r) is 
larger than 1.5. 

In fact, Eq. ^ is really valid for not too large a 
relativistic pressure in comparison to the thermal one 
(P re i /P th mu st be smaller than ~ 10, cf. Pelletier & 
Roland 1984) and we assume to be in this limit case so 



that we can also still neglect the thermal pressure of the 
plasma in comparison to the relativistic one. 

3.3.3. The high energy cut-off 

Since we take into accounts the cooling, a high energy 
cut-off in the particle distribution must necessarily appear 
at a Lorentz factor j c where heating and cooling balance 
(Webb et al., 1984). Since we suppose that particles cool 



via inverse Compton process, and assuming that the soft 
photon density is homogeneous in the shock region, the 
mean cooling time scale can be written: 

W = ?-^ = ?T L ^ (10) 



4^cr T W so ft 4Z S 7 c u{ 
where l s is the local soft compactness I 



oft^T 



AnRm„ 



being the electron mass,L so /t the soft radiation luminosity 
injected into the shock region and W so ft the corresponding 
local soft photon energy density. The maximum Lorentz 
factor 7 C achievable by the acceleration process is thus 
obtained by setting t acc = t coo i which gives, using Eqs. 
(I and ©: 
R u\ r — 1 



7c 



(11) 



4. Basic kinetic equations at shocks 

The shock region is depicted in Fig. [l] with the accelera- 
tion and pair creation region described above. The scale 

1 here s is the spectral index of the spatially integrated dis- 
tribution function of the particles 11(7) oc 7~ s . It is related to 
the spectral index q denned in Appendix by s = q — 2 



are not respected. We take the origin of the x-axis at the 
point where the second derivative of the flow velocity 
vanishes, that is: 



d 2 u 
dx 2 



0. 



(12) 



As previously said, we suppose the existence of a mag- 
netic field B, slightly perturbed by Alfven waves. The par- 
ticles are thus scattered by these waves trough pitch angle 
scattering (Jokipii 1976; Lacombe 1977) and can cross the 
shock front several times before escaping unless they are 
rapidly cooled by radiative processes. During these scat- 
terings, the particle gain energy trough the well known 
first order Fermi process. We also suppose the magnetic 
perturbations to have sufficiently small amplitudes so that 
we can treat the problem in quasilinear theory using the 
Fokker-Planck formalism. Besides we assume that, in each 
part of the shock front, the scattering is sufficient for the 
particle distribution to be nearly isotropic. With these 
different assumptions, and when first order Fermi pro- 
cess just as radiative losses and pair creation/annihilation 
are taken into account, the particles distribution function 
f(p, x) must verify the following equation : 



dt 



df 1 du df 1 db P i f 



dx 



3 dx^ dp 



d df 
— D— 

dx dx 



(13) 



p 2 dp 
+A(p)+B ± (p) + C±(p) 
The three first terms of the right member correspond to 
the first order process, the radiative losses (b being > 0, 
in the case of Compton cooling in the Thomson regime 
b = l/(jt coo i) where t coo i is given by Eq. |Io|) and the spa- 
tial diffusion respectively, A(p) is the injection of particles 
at 7 m m due to the pre-acceleration processes (cf. section 
|3.1.l| ), B±(p) is the pair creation rate and C±(p) the an- 
nihilation one. This equation is essentially the equation of 
the cosmic-ray transport originally given by Parker (1965), 
Skilling ( 1975| and references therein) except for the addi- 
tion of the radiative, pre-acceleration and pair processes. 



4-1- Evolution equation of the flow velocity profile 

As previously said in section ||, we supposed the shocked 
plasma pressure to be dominated by the pressure of the 
relativistic particles. The latter can thus be written as 
follows: 



P rel = J 4?r P 2 y /(P) x ) d P- 



Consequently, by multiplying Eq. ( |l3| ) by 4-7rp 3 c/3 and 
integrating on p we obtain the hydrodynamic equation 
linking pairs (through their pressure) and the flow velocity 
that is (in stationnary state): 



dP, 



el 



du d „ d 



dx ' 3 Te ^dx dx^ dx re ^ 



Q- +P± + Q+ (14) 



where: 
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corresponds to the pressure loss rate due to radiative losses 
(it is thus obviously negative), 

is the pair pressure creation rate due to the pre-accelerator 
processes, and 



4itp 6 c 



B±(p)dp- 



4irp 6 c 



C±{ P )dp 



3 ^ J r J 3 
is the pressure creation rate including the pair creation 
and annihilation processes. It is now possible to deduce 
from Eq. ( |l4]) the differential equation followed by the 
flow velocity. Indeed the momentum conservation equation 
gives: 



du 



du 



dR 



rel 



0. 



dx ) dx 

which, in stationnary state, is easily integrated: 

pu 2 + P rel = C (15) 

where C is a constant. Assuming the mass is dominated 

by protons (i.e. p = n p m p ) and thus is conserved, Eq. ( |l5| ) 

gives: 

-P re l = Piui(ux - u). 

To obtain this equation we have supposed f re l|_ = 
and «|_^, = u\. Then Eq. (|l4|) gives finally: 



Adu 
3 dx 



1 —00 
3 dx 



— 1 

dx 2 



d_ 

dx 



da 

dx 



n 



R~. 



where the reduced variables are defined as followed: 

~ u 



n 



X 

Q + 



and R~, 



P± + Q~ 



pml 



Li 



(16) 

(17) 
(18) 

(19) 



The parameter II is the ratio of the flux of energy trans- 
mitted to particles, through pair creation/annihilation, 
pre-acceleration and cooling processes, integrated on the 
"ID" volume -R 77 , to the flux of kinetic energy of the up- 
stream flow. As the acceleration and the cooling time are 
short with respect to the travel time i? 77 /c, we can as- 
sume that the sum of the 3 terms Q + , P±, and Q~ of Eq. 
( |l9| ) is equivalent to the instantaneous creation of pairs at 
the Lorentz factor 7 m m, that is: 

Q++P ± +Q- ~ 7 """ 3 Wec2 J 4np 2 B ± (p)dp 

2 

= ^ n± (20) 

where n± is the total (integrated over the particle Lorentz 
factor) pair creation rate. Then, LT can also be seen as 
the ratio of the pair luminosity (the pair power density 
integrated on the "ID" volume i? 77 ) to the kinetic energy 
flux of the flow. This parameter will play an important role 
in the evolution of the shock profile as we will see in the 
following. We can anticipate that a value of II of the order 
of the unity will be certainly unfavorable to the formation 
of the shock. It means that all the kinetic energy of the 
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Fig. 2. Plot of the 3 solutions of Eq. ( p2|) versus LT for 
different value of i? 77 . From thin to thick lines, i? 77 = 10, 
100, 10 3 . The flow velocity ttshock at the shock location, 
i.e. in x — is plotted in dashed line. For each value of 
i? 77 , there is a maximal value of LT, above which only one 
real solution exists. The disapperance of the ttshock branch 
means that the shock cannot exist anymore. 



upstream flow will be dissipated in radiative cooling and 
pair creation/annihilation processes. 



5. Solution without feedback 

In this section, we study the solutions of the previous equa- 
tions as a function of the parameter II. We do not consider 
the feedback of the hydrodynamics on the pair creation 
process, i.e. II is considered as a fixed free parameter. The 
complete problem, which must take into account the dif- 
ferent processes to find self consistent solutions, will be 
treated in the next section. 



5.1. A maximal pair creation rate 



A way to solve Eq. (|16|) is to integrate it between —00 and 
1/2- In this case, we can neglect the cooling during the inte- 
gration. Indeed, even if for x < —L\ particles do not inter- 
act with the shock, we have assumed that some processes 
apply and balance coolings so that particles are injected 
at a minimum Lorentz factor 7 m i„ above the resonance 
threshold 7 SCQ t. On the other hand, in the shock region 
(that is —L\ < x < L2), we have seen in section 3.3.3] that 
the coolings are negligible, in comparison to heatings, for 
particles with a Lorentz factor smaller than 7 C (the case 
of the majority of the particles). Besides, we can also ne- 
glect the annihilation since it occurs mainly at low energy, 
i.e. for particles with Lorentz factor 7 ~ 1 (cf. Coppi & 



Blandford 1990), whereas we assume 7 > 7 m i„ > 1. In 



these conditions, and with the assumption that the pair 
creation is homogeneous in the region 
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(meaning P± roughly constant) and null outside (meaning 
P±=0), the integration of Eq. ( |l6| ) gives: 



On 
c)x 



r(l-5) 



n 



R^ 



(21) u*{x) 



Equation ( pi| ) can be solved with the following boundary 
conditions: 

U{ — Ryy) — U*( — Ryy), 

where u* is the flow velocity solution of Eq. 

II = that is (Drury et al. |1982[ ): 

13 , [2 f 1 dx r 

tanh - / 

7 7 7.L D 



(Ell) when 
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and assuming Eq. (12). Using Eq. (J2l|), this last equation 
becomes: 
d 2 u 

— = P{u) = AQu 3 - Mu 2 + (39 + 42n)u 
ox 1 



-6n 4 



R. 







(22) 



This equation possesses in general 3 real solutions which 
obviously depend on II and i? 77 . They have been plotted 
in Fig. | as a function of II and for different values 
of Rjy, Of course, if II vanishes the three branches of 
solution converge respectively to the well-known results 
u=l/7, 4/7 and 1 corresponding to the values without 
pair creation. By continuity, the flow velocity at the shock 
location will follow the second branch noted Ugj^^ on 
the figure (plotted in dashed line). It appears that, for a 
given value of i? 77 , there exists a maximal value H max 
of II above which there is only one real solution which 
still verifies Eq. (p2|). The unphysical discontinuity of the 
flow velocity at the shock location for II = H max means 
simply that the shock cannot exist anymore. 

We have reported in Fig. |[ different shape of the 
polynomial P(u) and the corresponding flow velocity 
profiles obtained numerically by solving Eq. (|2l]) for 
different values of II (fixing i? 77 to 10). We clearly see 
the softening of the flow velocity profile when II increases 
meaning that the acceleration becomes less and less 
efficient. We have also plotted in Fig. 0, the variation 



of the compression ratio r (as defined in section 3.3.1) 
in function of the pair creation luminosity II and for 
different values of i? 77 . As expected, r converge to ~ 1 
when II increases. 



Since P(u) is a third degree polynomial, the transition 
between 3 to 1 real solution happens when the following 
conditions are satisfied: 
P(u) = 
P'(u) = 

The resolution of this system of equations gives thus a 
relation between H ma x and R^: 



63IL 



2/3 



R 



(23) 



We have plotted H ma x versus i? 77 in Fig. ||. We see that 
it is an increasing function of i? 77 , meaning that the 
larger the pair pressure creation region the larger the 
power we need to kill the shock. 

But the important conclusion of this part is that Umax 
is necessarily smaller than 3/14~ 0.20 (cf. Eq. (p3[)) mean- 
ing that at most 20% of the kinetic energy flux of the up- 
stream flow transformed in pairs is sufficient to suppress 
the shock discontinuity. 



0.05 0.1 0.15 

Pair luminosity II 



0.2 



Fig. 4. Variation of the compression ratio r in function of 
the pair luminosity II for different values of R 11 . i? 77 =10, 
100, 1000 for the solid, dotted and dot-dashed line respec- 
tively. 



0.25 



0.15 




0.05 - 



10 100 
Pair creation region size R,,/Li 



1000 



Fig. 5. Plot of the maximal pair luminosity Umax versus 
the size i? 77 = i? 77 /Li of the pair creation region. 



6. Stationary states 

6.1. The equations 

At this stage, we have studied the deformation of the 
flow profile due to the pair creation as a function of the 
parameter n. Conversely the velocity profile of the flow 
will control the distribution function of the accelerated 
particles, and hence the number of particles able to trigger 
the pair creation process (i.e. particles with 7 > 7th)- 
This back reaction will thus select the control parameter 
II defining the possible stationary states for a given set of 
external parameters. 

From the previous section we have seen that the 
compression ratio r decreases when the pair creation 
rate raises (cf. Fig. ||) but reversely if r decreases, the 
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spectral index will increase and the number of high 
energy particles will decrease (cf. Eqs. (||) and (pi])), the 
final effect being a decrease of the pair creation rate. 
Consequently, we indeed expect the system to reach, in 
some conditions, stationary states where hydrodynamics 
and pair creation effects balance. 

The problem is now to calculate the pair production 
rate as a function of the flow profile parameters. For this 
goal, we first approximate the downstream particle dis- 
tribution function just after the shock with the following 
form : 

2. 

Re- 
appears to be a 



712(7, x = ) = n ol s exp 



(24) 



From the remarks of section 3.3 



it 



relatively good approximation in the vicinity of shocks 
where first order Fermi process and Compton coolings 
occur. The assumption of an exponential high energy 
cut-off appears also to be reasonnable (cf. for cxcmplc 
Webb et al. |98l. 



The particle energy distribution is then completely de- 
termined by no, s, and j c . The latter are function only 
of the compression ratio r and the external and geomet- 
rical parameters i?, l s and u\/c. We can thus express 
12(7, x = + ) as a function only of these external param- 
eters and the two dimensionless variables r and 

~ TTlp 

■no (25) 



"0 



Pi 



The compression ratio r is a function of II and i? 77 
through the numerical integration of Eq. ((2^). We can 
express also no as a function of these two quantities in the 
following way. For simplicity, we will neglect, in the ac- 
celeration processes, particles coming from the upstream 
flow and consider the plasma to be dominated by created 
pairs. Consequently, assuming that the pair creation rate 
h± is constant within the pair creation region, the total 
downstream particle density n to t after the shock can be 
written: 

R 

n 2 {i,x = + )d 7 ~ n± — ^. (26) 



n to t 



With the help of Eqs. <M) and fl2 



77 

we thus find: 



in,. 



Pi 



-n±- 



R 11 (s — 1) 



U 2 



where we 



7 exp 



Imin 

have 

2 



use 
7 1_s 

/ mm 

' s-r 



(27) 

the approximation 
reasonable for s > 1, 



which is always our case (cf. section 3.3) 



As we assume that the produced pairs are instanta- 
neously accelerated to the Lorentz factor jmim the pair 
luminosity II can be approximated by (cf. Eq. (|T^) and 

©): 

— . (28) 



n 



Consequently, Eq. (E71) can be re-written: 



77,0 — II— jrity-y- 

c z 



1) 



7 2 " s 

' ru i n 



(29) 



which links no to Lf and R 



77" 



To close our systems^of equations, we now need to link 
the parameters II and R 11 to no and r using two other 
independent equations. This could be done by comput- 
ing directly the pair creation rate, as follows. The parti- 
cles accelerated in the shock cool via IC process on the 
surrounding external soft photons. The photon spectrum 
emitted by a single particle of Lorentz factor 7 per unit of 
time and dimensionless energy e can be approximated by: 
dN ph 



dtde 



c l s 1 da c l s 4 2 

o T = ~5 — °( e ~ Cs > 

R e s ar de Re s 3 



where we have approximated the differential Compton dif- 
fusion cross section by a Dirac distribution, i.e —r- = 

de 

a T 5(e~ -7 e s ). 

The photon spectrum emitted by a particle of initial 
Lorentz factor 7 cooling down to 7 = 1, is then given 
by: 



dN ph 



de 



dN 



ph 



dt 



dtde 



dt 



where the cooling rate — — 



. . „ ~ — o 7 / 1 ' -i^uu 

ay 7 

(|l0|). The integration of Eq. (|g) then gives 



(30) 



, with tcooi given by Eq. 



dN, 



ph 



de 



I6e< 



H{\l 2 e s )] 



where H is the Heaviside step function. It is now possi- 
ble to compute the total photon spectrum emitted by the 
distribution of the particles crossing the shock per unit of 
time and dimensionless energy e: 

dN ph 



dNX 



dtde 



(e) 



.Alport (7) 



dc 



where N part (j) = ^(7,2; = + )u2irR 2 is the flux of par- 
ticle crossing the shock. Consequently: 



dtde 



(0 



nohin c 2 Rc 



x exp 



u\ a T 



16e s 



-£±2 1th 



7c 



(31) 



where we define the kinetic compactness Ikin = 
(piv$irR 2 )a T , , 

5 . We can then deduce the pair creation rate 

Rm e c i 

due to the photon-photon annihilation on a length scale 



1 



n± 



ttR 2 R., 



dtde 



(e) (l 



L 77 J e=l 

where T 77 (e) is given by Eq. ( 
relation between II, n and r: 



- e- T ^j de (32) 
We finally obtain a new 



n 



— 2 2 — 
Tmin^O ^ ^th 

4r ui S — 
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exp 

7c 

Finally, i? 77 is given by Eq. 
L.- 3 



1 _ e -^iW),/ ( 
) where: 



--2) 
1 



(33) 



(34) 



dN to h l 

p -(e). Together with Eqs. @, 



^ ""^ = 7T./?-f ,/,,/, 

(E9) and (B3J) , this forms a system of four equations 
of four unknowns II, i? 77 , hq and r which has to be 
verified at the equilibrium i.e. when pair pressure and 
hydrodynamics effects balance. 



6.2. The solutions 

The above system depends on six different parameters: 
7mm j the minimal Lorentz factor for the particles to 
be accelerated in the shock, e s , the soft photon energy 
(in unit of m e c 2 ), u\, the upstream flow velocity, R 
the transverse size of the shock (in unit of the diffusion 
length Li), l s , the soft compactness and Ikm the kinetic 
compactness. 

The system is solvable in only some part of the 
parameter space. It always possesses two solutions: a 
"high pair density" one (large n , large II, small r) and a 
"low pair density" one (small no, small II, large r). The 
latter connects to the trivial solution of the problem i.e. 
II = 0. However, only the high pair density states are 
compatible with the neglecting of particles coming from 
the upstream flow (cf. Eq. ^6|). Thus we will only focus, 
in the following, on this branch of solutions. 



tion of r (cf. Eq. |Tl|) , which decreases if r decreases (cf 
Fig. |b). 

— Variation of "f m in' since we suppose that the pairs are 
instantaneously accelerated to the Lorentz factor 7 mm , 
the larger 7 mm and the larger the pair pressure i.e. the 
smaller the compression ratio (cf Fig. |^c). We recall 
however that for too small values of 7 mm (i.e. jmin — 
1) the annihilation of pairs is not negligible anymore 
and our model is no more valid. 

— Variation of l s : the larger the soft compactness, the 
shorter the cooling time scale (cf. Eq. (|ic|)) and thus 
the smaller j c . An increase of l s thus disfavors the pair 
creation process. The compression ratio has then to 
increase to keep the pair-hydrodynamics balance. For 
small l s however, j c is very large (3> Jth) an( i its value 
becomes immaterial. Consequently, II and r, become 
independent of l s , as shown in Fig. |^d. 

— Variation of u±: when u\ increases, t acc decreases and 
then 7 C increases. The compression ratio r decreases 
by compensation. On the other hand, for low values of 
u\ (i.e. u\ < 0.1c), 7 C becomes too small, decreasing 
the pair creation rate and r increases rapidly. 

— Variation of I kin- f° r high values of I kin, the pair den- 
sity must be also high to efficiently modify the hydro- 
dynamical profile. The pair creation process is thus sat- 
urated meaning that the pair creation rate grows lin- 
early with Ikini i-e. II, r and j c keep constant. However, 
for low values of Ihini the pair density is so low that 
the pair creation optical depth becomes smaller than 
1. In compensation r must increase, so n decreases 
accordingly. 



For clarity, and in order to better understand the influ- 
ence of each parameter on the solutions of our problem, we 
have plotted in Fig. ^, the compression ratio r, the cut-off 
Lorentz factor 7 C and the pair luminosity n versus each 
one of these parameters, the other being fixed to values 
indicated in the figure caption. Some comments may be 
done on these different figures: 

— Variation of R: an increase of R will produce an in- 
crease of i? 77 (cf. Eq. ||) and thus an increase of n 
and a decrease of r. For R ^> L a , i? 77 converges to 
L D so that n and r become constant (cf Fig. ||a) . The 
cut-off Lorentz factor 7 C increases with R as expected 
from Eq. ([ll]), because a constant l s implies a decreas- 
ing photon density. For low values of R, j c become 
smaller than j t h and the decrease of the pair creation 
rate must be compensated by a hardening of the dis- 
tribution (r increases and n decreases). 

— Variation of e s : when e s increases, ^th decreases, favor- 
ing the pair creation process (n increases). It is com- 
pensated by a steepening of the particle distribution 
(i.e. r decreases). Concerning 7 C , it is a simple func- 



7. Discussion 

In this paper, we have seen that the launch of the pair 
creation process by particles accelerated by a shock em- 
bedded in a dense soft photon field, could disrupt, through 
the increase of the associated pair pressure, the plasma 
flow and eventually, for too high pressure, to smooth com- 
pletely the shock. By including the feedbacks of the hy- 
drodynamic of the flow on the pair creation process, the 
system can reach, in some conditions, stationary states. 
Such processes may be at the origin of the high energy 
emission observed in compact objects where acceleration 
processes in dense photon field are expected to occur. We 
give here numerical estimates of the physical parameters 
relevant to this class of objects. 

7.1. High energy spectra 

As seen in the previous section, for a given set of parame- 
ters, the system may reach stationary states character- 
ized by a compression ratio r. The spectral index and 
cut-off of the particle distribution function just after the 
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0.0 0.2 0.4 0.6 0.8 1.0 10 2 10 3 10 4 10 5 



Upstream flow velocity u,/c Kinetic compactness l kln 

Fig. 6. Plots of r, II and 7 C versus each parameter of the model: a) the reduced shock transverse size R, b) the reduced 
soft photon energy e s , c) the minimum Lorentz factor j m i n , d) the soft compactness l s e) the upstream flow velocity 
Ui/c and f) the flow kinetic compactness the other parameters being fixed to R — 10 6 , e s = 2 x 10~ 5 , j m in = 10, 
ui/c = 0.1, l s = 100 and l km = 10 4 . 



shock n 2 {j,0 + ) are then given by Eq. (||) and (|T^) re- 
spectively. Since we have assumed that these particles are 
cooled via IC process (where the cooling length scale de- 
pends on the particle Lorentz factor), the emitted energy 



spectrum is characterized by a cut-off power law shape 
E 



Fe oc E a exp 



E r 



(cf. Eq. (pll)) where a and 
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Fig. 7. Contour plots of the spectral index a (solid lines) and the high energy cut-off E c in keV (dashed lines) of the 
emitted spectrum in the (l s , Ikin) space. The soft photon energy is equal to 10 eV and 600 eV for the left en right 
plots respectively and R = 10 4 and 10 6 for the lower and upper plots respectively. The other parameters have been 
fixed to 7mm = 10 and U\/c — 0.1. In the grey region there is no pair dominated solution in stationary state . 



(35) 
(36) 



E c are simple functions of s and 7 C : 

s 

a = 2 

We have reported in Fig f^the contour plot of a(l s , Ikin) 
and E c (l s , Ikin) for two values of e s (10 and 600 eV which 
are representative of the typical values of soft photons 
emitted by an accretion disk around a supermassive and 
stellar mass black hole respectively) and R (10 4 and 10 6 ). 
The other parameters have been fixed to J m in = 10 and 
U\/c = 0.1. The contours of a and E c keep roughly the 
same shape but cover a different region of the parameter 
space for different parameter sets. In each figure, we can 
see that: 



The spectral index does not strongly vary between the 
harder and the softer spectra (Aa ~ 0.3). It reaches 
an asymptotic plateau for low l s and high lki n ■ In these 
conditions, both pair density and high energy cut-off 
are large so that the pair creation process is saturated. 



Then it can be shown that Eq. ( p2| 
tively simple relation between a, 7. 
1— a 



reduces to a rela- 
mm and e S) namely: 



1 

2a 



3 
4^ 



Concerning E c , following Eqs. ( [TT| ) and (|3^), it is in- 
versely proportional to I 2 . 

— The harder spectra are obtained for large values of l s 
or small values of lkin ■ In both case, the pair efficiency 
decreases either due to a small value of E c or a low 
particle density, that is a low pair production optical 
depth. To keep the equilibrium between the pair cre- 
ation effects and the hydrodynamic of the flow, the 
system has to reach harder spectra to compensate this 
decrease of the pair creation rate. 

— For small l s (respectively large Ikin), ct is independent 
of l s (respectively Ikin)- This is a direct consequence 
of the behavior of the compression ratio at small l s 
or large Ikin as explained in section S.2 
that in the case e 



600 eV and R 



We can note 
10 6 shown in 



Fig. (7] the condition j c < m e c 2 /e s , needed in Thomson 
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regime (cf. section 2.5), may not be verified in the low 
part of the parameter space. However, since the pair 
creation process is mainly triggered by particles with 
Lorentz factor of the order of jth (<C 7 C ), we do not 
expect the Klein-Nishina effects to strongly modify our 
results. 

— For to large l s or to small the system cannot reach 
sufficient hard states to keep in equilibrium and no 
high pair density stationary state can exist any more. 
The system can only be in the trivial pair free state 
(i.e. r = 7 and II = 0). 

We interpret the differences between the four plots of Fig. 
[?] as follows. The increase of e s favors the pair creation pro- 
cess. Thus, for given values of Z s and Ihini the spectral in- 
dex is larger. The high energy cut-off E c increases mainly 
because of its dependence on e s (cf. Eq. |36| ), j c keeping 
roughly constant (cf . Eq. |ll|) . Concerning the variation of 
R, as expected, a larger value corresponds to a more effi- 
cient acceleration (cf. Eq. [ll]) and favors the pair creation 
process which allows a larger space parameter region for 
high pair density solutions. 

7.2. Annihilation line 

The presence of pairs should give a signature as an 
annihilation feature at ~ 511 keV. We will show that 
this feature is not expected to be strong in our model. 
Here, we have supposed the existence in the shock region 
of pre-accelerating processes bringing leptons to the 
sufficient energy (i.e. 7 > j m i n with j m i n of the order 
of a few, cf section |3.1.l| ) for resonant scattering off 
magnetic disturbances. Since the annihilation process 
occurs mainly at low energy, i.e. for particles with Lorentz 
factor 7 ~ 1 (cf. Coppi & Blandford 1990), it occurs 



the luminosity of the annihilation radiation to be smaller 
than few percent of the total high energy radiation, which 
is quite compatible with the non observation of strong 
annihilation lines in this class of Seyfert galaxies as shown 
by the best upper limit observed in Seyfert galaxies with 



the OSSE satellite (Johnson et al. 1997) 



mainly far downstream, where the pairs created in the 
shock can cool down. The annihilation line luminosity 
is thus at most equal to the pair rest mass luminosity, 
which is itself smaller by a factor j m in than the pair 
creation luminosity (Eq. ^8|). As shown in section |5.l| , 
the pair luminosity n is itself limited and is necessarily 
smaller than ^20% of the X-ray/7-ray luminosity, i.e. 
n < Umax — 0.2 (assuming that the total kinetic energy 
of the upstream flow is transformed in radiation). Besides, 
for n ~ 0.2, the compression factor is very small, of the 
order of unity, resulting in a very steep X-ray spectra. 
When hydrodynamics feedback is taken into account, the 
pair luminosity may be well below this theoretical limit 
of 20%. An X-ray spectral photon index ~ 2 (as those 
generally seen in Seyfert galaxies) requires a compression 
ratio r ~ 3-4. Such values of r require values of n smaller 
than ~ 10% (cf. Fig|). 

Assuming a steady state, pairs annihilate at the same 
rate as they are produced. H/j m i n gives then an upper 
limit of the annihilation line luminosity. We thus expect 



7.3. Variability 

For some values of external parameters, as suggested by 
Fig. ^, no high pair density solution can exist in station- 
ary states. So only the pair free solution (i.e. with n = 0) 
exists. The system is not expected to be variable with a 
constant set of parameters and variability can only occur 
with a variation of one of them. An interesting possibility 
would be to consider a possible feedback of the relativis- 
tic plasma to the soft compactness. In the reillumination 
models for instance (Collin, 1991; Henri & Petrucci, 1997), 
the soft photons are produced by the reprocessing of the 
primary X-ray emission. An increase of the pair plasma 
density will increase the X-ray illumination and thus the 
soft compactness. Fig. fj] shows that in some cases, the 
change of l s make the system switch to pair free solution 
which will stop the reillumination and bring the system 
back to pair rich solutions. Limiting cycles could thus oc- 
cur. We intend to further investigate this possible effect 
under astrophysically relevant conditions. 



8. Conclusion 

In the present paper, we have studied the effect of pair 
creation, via high energy photon-photon interaction, on 
a shock structure, where the high energy photons are 
produced via IC by the particles accelerated by the shock 
itself. The problem is highly nonlinear since pairs can 
modify the shock profile through their pressure and, 
mutually, a change of the shock hydrodynamics can 
decrease or increase the pair production rate. 

We have shown that for a given size of the pair 
creation region, it exist a maximal value of the pair 
creation rate above which the shock cannot exist any- 
more. When the hydrodynamical feedbacks on the pair 
creation process are neglected, a pair power of at most 
20% of the upstream kinetic power is sufficient to kill 
the shock. This constraint can fall to few percents in 
stationary states where pair creation and hydrody- 
namical effects balance. We thus do not expect the 
presence of strong annihilation lines. We also obtain spec- 
tral parameters in rough agreement with the observations. 

We suggest also a possible variability mechanism 
if the soft photon compactness depends itself on the 
pair density of the hot plasma, such as expected in 
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reillumination models. 

In the model presented here, the cooling of particles is 
due to the IC process on external soft photons. However, 
particles may also cool on soft photon they produce by 
synchrotron process when spiraling around the magnetic 
field lines (the so-called synchrotron-self-compton process, 
SSC). In this case the cooling will also depend on the parti- 
cle distribution function. We may expect that the addition 
of the SSC process would allow to obtain stationary states 
with harder spectra than those we obtained here, since the 
additional synchrotron cooling would be compensate by a 
stronger acceleration, i.e. a larger compression ratio r. The 
detailed study of this problem is left to future work. 
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Appendix A: Particle distribution function with 
pair creation 

Neglecting the annihilation rate which is only important 
at low energy, and asssuming that pairs are produced at 
a constant rate Q in the region — i? 77 < x < +i? 77 and 
inst antane ously accelerated at p m in = IminmeC 2 (cf. sec- 
tion ft.l.lj ), Eq. (13) can be rewritten (in the stationnary 
state) as follows: 

1 du df 1 dbp 4 f 
3 dx^ dp p 2 dp 
Q5{p - p m in)H(R 17 - \x\) (A.l) 
where H(x) the Heaviside function defined by 
H(x) = 1 for x > 

= everywh ere els e. 
As shown in section 3.2.2] , the shock thickness (which is of 
the order of the diffusion length) is largely smaller than 
any of the characteristic leng th o f the problem. Conse- 
quently, when integrating Eq. (A.l) between — oo and + , 
we can still approximate the velocity gradient by a Dirac 
du 

-Q- = (u2 — u\)S(x). As we can also neglect the coolings 

for x < for particles with Lorentz factor 7 < -f c (cf. 
section |3.1.2 ). The result of the integration gives: 

r ( r ~ M d/ Q + t , QRry x , s 
/o+ + I ^7— ) = h H — - — o{p~p min ) 



3r 



dp 



Ul 

Ddf 0+ 



(A.2) 



ui dx 

where / + and f± are the distribution function in + and 
in — 00 respectively. In the downstream flow only the pair 
creation process will modify the distribution particles, and 
the main result will be an increase of the total number of 
particles on a scale of the order of i? 77 . We can thus make 
the following approximation: 
D df + D u 2 f + L 2 f + 



ui dx U2 ui Rf\ 
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Since we have i? 77 ^> L2, this term is negligible in 



Eq. 



), and its integration gives: 



fo+ (p) = OP 



MP') 



QR-y 
Ui 



'-S(p' -p m in) 



P 



,q-l 



dp' 



With Q = 0, one finds again the solution without pair 

creation. For p > p m i n = jminirieC, the integral is almost 

constant. We thus find that the distribution function is a 

3r 

power law of spectral index q = in any case. 

r — 1 



